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Dynamics of regio-regulated Poly(4-methylthiazole-2,5-diyl) [HH-P4MeTz] was investigated by 
solid-state 1 H, 2 D, 13 C NMR spectroscopies, and differential scanning calorimetry(DSC) measure- 
ments. DSC, 2 D quadrupolar echo NMR, 13 C cross-polarization and magic-angle spinning(CPMAS) 
NMR, and 2D spin-echo(2DSE) CPMAS NMR spectroscopy suggest existence of a quasi-ordered 
phase in which backbone twists take place with weakened 7r-stackings. Two-dimensional exchange 
2 D NMR(2DEX) detected slow dynamics with a rate of an order of 10 2 Hz for the CD3 group in 
d3-HH-P4MeTz at 288K. The frequency dependence of proton longitudinal relaxation rate at 288K 
shows a tj _1//2 dependence, which is due to the one-dimensional diffusion-like motion of backbone 
conformational modulation waves. The diffusion rate was estimated as 3±2 GHz, which was approx- 
imately 10 7 times larger than that estimated by 2DEX NMR measurements. These results suggest 
that there exists anomalous dispersion of modulation waves in HH-P4MeTz. The one-dimensional 
group velocity of the wave packet is responsible for the behavior of proton longitudinal relaxation 
time. On the other hand, the 2DEX NMR is sensitive to phase velocity of the nutation of methyl 
groups that is associated with backbone twists. From proton Ti and T2 measurements, the acti- 
vation energy was estimated as 2.9 and 3.4 kcal/mol, respectively. These were in agreement with 
3.0 kcal/mol determined by M0ller-Plesset(MP2) molecular orbital(MO) calculation. We also per- 
formed chemical shielding calculation of the methyl-carbon in order to understand chemical shift 
tensor behavior, leading to the fact that a quasi-ordered phase coexist with the crystalline phase. 

PACS numbers: 82.35.Cd,82.56.-b,76.60.-k 



I. INTRODUCTION 

In the discussion of structure-property correlation in 
polymeric condensed materials such as semicrystalline 
polymers and polymeric liquid crystals, one must of- 
ten take into account of not only duality of crys- 
tal/amorphous but also elemental excitations of poly- 
meric condensed systems or their relaxation from ele- 
mentally excited states to a ground state. For instance, 
exciton generation and anihilation in 7r-conjugated poly- 
mers and dynamic disorder in polymer crystals are typ- 
ical elemental excitations found in polymeric systems. 
"Conformon" is one of elemental excitations concerning 
molecular structured It is known that some of semi- 
crystalline polymers form mesophases that are called a 
quasi-ordered phased. The phase is defined as a partially 
disordered crystalline phase, of which an order parameter 
is larger than that of amorphous but smaller than that 
of perfect crystal. One may consider this phase as liquid 
crystal(LC) with high viscosity. 

Since the quasi-ordered phase of 7r-conjugated poly- 
mers is closely related to conjugation length of n elec- 
trons in a crystal, it is quite important to investigate 



both static and dynamic structure. In spite of the im- 
portance, information about the dynamics of the phase 
is still poor, mainly because there are quite limited ex- 
perimental methodologies to approach the problem. As 
a theoretical approach, van der Horst and co-workersi 
have proposed that in crystalline conjugated polymers, 
the presence of static disorder brings about a spread 
of the distance between two chains constituting the de- 
crease of wave-function overlap and the exciton induced 
by thermally excited intermolecular phonon modes (dy- 
namic disorder) will reduce the quantum mechanical co- 
herence. Hence, the conception of a single polymer chain 
put into a dielectric medium can interpret the features of 
conjugated polymers. In this case, the polymers have a 
character of one-dimensional substances because electri- 
cal conduction is dominated by the free electron having 
only degrees of freedom which can move in one dimen- 
sion. 

Heat capacity, quasielastic neutron scattering, and nu- 
clear magnetic resonance (NMR) experiments have inten- 
sively been performed to explore collective dynamics of 
quasi-ordered structure including charge and spin den- 
sity waves^, domain wall dynamics^ and discommensu- 
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ration in incommensurate systems 7 ". So far quasi-ordered 
phases in polymeric systems are, however, not yet well 
established. 

NMR spectroscopy is one of promissing methods to 
obtain information about molecular dynamics since one 
can perform nondestructive experiments. Particular at- 
tention has been devoted to the phase transiton of tt- 
conjugated small molecules which have a class of quasi- 
ordered phases, the incommensurate(IC) phase (such 
as biphenylSi^iill, p-terphenylii* 1 ^, bis(4-chlorophenyl)Ai 
and dichlorobipenylsulphone 14 , etc.). In these phases, 
disordered crystals show the long range breaking of the 
lattice periodicity. Therefore the wavelength of the mod- 
ulation is not an integral multiple of the unit cell, which 
means that translational symmetry is lost. In most in- 
commensurate crystals, power-law temperature and fre- 
quency dependences on proton longitudinal relaxation 
time were observed. Furthermore, the relaxation rate 
rises dramatically around the transition temperature, 
T c . This is due to the fact that the spin relaxation is 
controlled by low-frequency excitations of the collective 
phase modes (phasons) corresponding to the critical slow- 
ing down of molecular motion accompanying the phase 
transition. 

Dynamics of modulation waves in IC solids is also in- 
vestigated by various NMR experiments. In Rb2ZnCl4 
(and its analogue) 7 i 15 i 16 i 17 i 18 i 19 i 20 , bis(4-chlorophenyl) 
sulfone 13 , and barium sodium niobateSLSLS, it is known 
that the structural modulation wave with respect to 
atomic displacements from equiblium positions produces 
a discommensuration in the higher temperature regime 
of incommensurate-commensurate transition tempera- 
ture. As far as 7r-conjugated polymers are concerned, 
there are, however, quite limited publications concern- 
ing quasi-ordered phases. In irans-polyacetylene one- 
dimensional diffusion of electron spin along the chain 
has been observed 24,25 - 26 ' 27 ' 28,29 . This result could be 
described by using a well-known soliton (domain-wall) 
model which indicates that in the presence of defects or 
impurities the mobile spins are trapped and act as the 
domain wall. 

Poly (thiophene-2 , 5-diyl) [PTh] , poly (pyrrole- 2 ,5- 

diyl) [PPyrr] , and their derivatives (e.g. poly ( 3 - 
alkylthiophcnc -2,5- diyl [P3RTh] ), which are made up 
of recurring five-membered rings, have been extensively 
in";;tigatc d 30 ' 31 ' 32 ' 33 ' 34 ' 35 ' 36 ' 37 ' 38 ' 39 ' 40 ' 41 ' 42 ' 43 ' 44 ' 45 ' 46 ' 47 ' 48 ' 49 
because they are thought to take coplanar structures and 
to eventually form highly extended 7r-electron conjugated 
network^ 5 ^ 3 - 5 ^ owing to their less sterically hindered 
structure compared with those of poly(arylene)s with 
six-mcmbcrcd rings such as poly(|)-phenylcnc)[PPP] 55 . 

PTh and PPyrr are constituted of " electron-excessive" 
heterocyclic units and are susceptible to chemical or elec- 
trochemical oxidation, i.e. p-doping. On the other hand, 
it has been recently reported that 7r-electron conjugated 
polyheterocycles with five-membered rings contain- 
ing " electron- withdrawing" imine nitrogen(s), poly (4- 
alkylthiazole-2,5-diyl)[P4RTz], were synthesized with 



dehalogenative organometallic poly condensatio n 44 ! 45 ' 56 
and have interesting redox and optical properties^ 7 -. For 
solutions of P4MeTz and P3MeTh, there are no signif- 
icant differences in absorption maximum (A max ) of UV- 
visiblc optical absorption spectr Amax = 425nm for 
the trifluoroacetic acid(TFA) solution of regio-regulated 
head-to-head(HH)-P4MeTz and 420 nm for the formic 
acid solution of rand-P4MeTz. For the solvent-cast films 
of HH-P4MeTz, on the other hand, the value of A max 
(498nm) shows significant red shift compared to that 
for the film of rand-PiMeTz (420nm). Further, HH- 
P4MeTz gives relatively large third-order non-linear op- 
tical susceptibility, x^ 3 \ °f 2.5x 10 _11 esu. The value of 
%( 3 ) is eight times larger than that for ra?id-P4MeTz 
(0.3x 10 _11 esu). These regio-regularity dependence of 
the optical properties can be observed only for P4MeTz. 
There are no similar observations on P3RTh. 

The x-ray diffraction study indicates that HH-P4MeTz 
takes a face-to-face 7r-stacking structure while head-to- 
tail(HT)-P3MeTh takes the staggered tt- stacking 5 ^. 
One of the interesting structural features of HH-P4MeTz 
is an alternative layered structure; one is a layer consti- 
tuted by highly densed methyl groups with two dimen- 
sional manner and the other is built by the face-to-face 
7r-stacking (Fig^l. 

In this article, we have synthesized natual abun- 
dant and methyl-deuterated (perdeuterated) versions of 
HH-P4MeTz, and have performed differential scanning 
calorimetry(DSC), variable temperature 2 D quadrupolar 
echo NMR, two-dimensional exchange (2DEX) 2 D NMR, 
variable frequency and temperature proton longitudinal 
relaxation time (Ti), and variable temperature proton 
transverse relaxation time(T2) measurements, in order 
to explore the dynamics of the quasi-ordered structure 
of the polymer. From these experiments, we shall dis- 
cuss wave length dispersion of modulational waves in the 
quasi-ordered phase of the polymer and confirmed that 
there exist the strongly anisotropic diffusion and anoma- 
lous dispersion of the structural modulation waves. 
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II. EXPERIMENTAL 



Materials 



HH-P4MeTz (M W =320G; light scattering) was pre- 
pared by the previously reported schemes involving 
dehalogenative organometallic polycondensationi 5 .. A 
methyl-deuterated version of HH-P4MeTz was synthe- 
sized from a scheme initiated from mono-bromination of 
deuterated aceton. The detailed scheme for the synthesis 
will be described elsewhere. For all the measurements, 
we used powder samples of HH-P4MeTz, which are re- 
crystallized by hexanuoroisopropanol(HFIP). 
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B. Differential scanning calorimetry measurements 

DSC thermograms of natural abundant and perdeuter- 
ated versions of HH-P4MeTz were recorded on a Seiko 
DSC 220 system connected to a SSC5300 workstation. 
The samples were first heated from 173K to 423 ~ 473K. 
The heating rate were 20 K/min. (HH-P4MeTz) and 
5, 20 K/min. (d 3 -HH-P4MeTz). After the first heat- 
ing scans, the samples were rapidly quenched by liquid 
nitrogen and then they were again heated at the same 
rate. 



dimensional techniques. The 2DEX NMR method partic- 
ularly provides the valuable data when trying to observe 
slow dynamics in solids 20 - 23 A spatial motion of 
certain atom (or molecule) such as diffusion or reorienta- 
tion is detected in the form of frequency change during 
the mixing time, t m j X . Moving nuclei experience a change 
of resonace frequency between at the beginning and at 
the end of t m ; x and then create cross peaks, while static 
nuclei without changing the frequency (or moving nuclei 
with a frequency of much higher than l/t m i x ) produce 
the diagonal signals in the 2D spectrum. 



Solid-state NMR measurements 



3. 



3 C NMR spectroscopy 



1. 1 H solid state NMR spectroscopy 

We carried out variable frequency proton longitudinal 
relaxation time (Ti) measurements in order to investi- 
gate the spectral density function of fluctuation of the 
local field in HH-P4MeTz. The proton longitudinal re- 
laxation times were measured using following three meth- 
ods; the saturation-recovery method(Tin) for three Lar- 
mor frequencies of 25(0.59 Tesla), 270(6. 34 Tesla), and 
400MHz(9.4 Tesla), the longitudinal relaxation time in 
the rotating frame of a radio frequency field(Ti p )£i, and 
in the dipolar order built up by the Jeener-Broekaert 
method(TiD)^ were measured at a resonance frequency 
of 400MHz. For T% p measurements on a wide range of 
radio-frequency (rf) field, we made a home-built probe 
with a solenoid type microcoil with a radius of 1.2mm 
and a length of 2.5mm. The intensity of rf field was in 
the range from about 43kHz to 1.37MHz. 

The proton transverse relaxation time(T2) measure- 
ments were performed at a resonance frequency of 
400MHz. The two-pulse Hahn echo (TPHE) and 
the Carr-Purcell Meiboom-Gill(CPMG) methods can be 
used to observe diffusion- like motions under the condition 
of local magnetic field gradient in a crystal^ 3 ". The experi- 
ments were performed as a function of echo time. If there 
is no diffusion, the apparent T2 by CPMG shows a fixed 
value independent on the echo time. On the contrary, 
when diffusive motion exists the value of T2 is changed 
with a variation of echo time. 



13 C cross-polarization and magic- angle sample spin- 
ning(CPMAS) NMR measurements were performed on 
a JEOL GSX-270 FT NMR spectrometer (6.34 Tesla). 
Solid-state two dimensional spin-ccho(2DSE) 13 C CP- 
MAS NMR measurements were performed on a Varian 
Unity400 FT NMR spectrometer (9.4 Tesla) equipped 
with a double resonance tunable MAS probe. The rf field 
intensities during the Hartmann-Hahn cross-polarization 
were set at 50kHz for the 1 H and 13 C channels. A contact 
time for the cross-polarization was 2.0ms. The intensity 
of CW proton decoupling was 64kHz. The recycle de- 
lay for all the experiments was 5s. The MAS speed was 
monitored and controlled by a personal computer with 
optical fibers. Principal components of chemical shift 
tensors for methyl carbons were determined with 2DSE 
spectroscopy, which is originally developed by Kolbert 
et al£L and applied to determination of CST compo- 
nents with small CSA by Asakawa et al~. A pulse se- 
quence of (preparation) —txjl — tt — t\/2— acq^) was 
used for 2DSE measurements. The detailed setting up 
for the 2DSE experiments were the following; the MAS 
speed was set at cj i /2tt = 1025Hz ± 5Hz, 128 t\ values 
were collected, with 512 acquisitions per t\ value, and the 
time increment in the t\ dimension was 122/is (the incre- 
ment for each t\j1 duration was set at the rotor cycle of 
1/16). 



III. THEORETICAL 



2. 2 D NMR spectroscopy 

We used a standard quadrupolar echo pulse sequence 
for 2 D NMR measurements. The rf-field intensity was set 
at 64kHz. The recycle delay for all the experiments was 
30s. Signal averagings were performed with 32-256 scans. 
The echo time for quadrupolar echo measurements was 
set at 0.020ms. 

For investigation of dynamic processes such as chemi- 
cal exchange, cross relaxation, spin diffusion, and atomic 
hopping motion, the two-dimensional exchange spec- 
troscopy has a number of advantages over the one- 



A. Spectrum Simulation 

A NMR spectral simulation program for 2DSE was 
written by C language and were performed on an IBM- 
AT compatible personal computer with using a GNU C 
compiler. Powder averagings were performed with ran- 
dom orientations with respect to the external magnetic 
field. The best fit simulation was picked out by monitor- 
ing the value of e, which is defined as 

_ -I (Em En Jmnlmn) 2 

V V J 2 V V P ' 

Z-^va Z-~tn ^ van L^dva Z-~in ran 
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where J mn and 7 mn are the experimental and simulated 
signal intensity for the (m,n)th spinning sideband. To 
describe a chemical shielding tensor, we used the span(fi) 
and skew(K)^ 9 . as well as the S and r\. 

B. Chemical Shielding Calculation 

13 C chemical shielding tensors of HH-P4MeTz were 
calculated by the ab initio self consistent field(SCF) cou- 
pled Hartree-Fock method with gauge invariant atomic 
orbitals(SCF-GIAO)22i 7 i and the second-order M0ller- 
Plesset GIAO(MP2-GIAO)S. Tail-tail(TT)-bi(4MeTz) 
were employed as model compounds for all the calcula- 
tions and were optimized by the MP2 method with the 
6-31G(d) basis set, and the SCF-GIAO and MP2-GIAO 
shielding calculations were carried out with 6-3 1G basis 
set. All the ab initio chemical shielding calculations were 
performed with Gaussian 98(Rev.A7) program package 7 ^ 
run on a Cray C916/12256 super computer at the Com- 
puter Center, Tokyo Institute of Technology, and SGI 
Origin2000 or Fujitsu VPP2800 at Institute for Molecu- 
lar Science in Okazaki, Japan. 

IV. RESULTS AND DISCUSSION 

A. DSC measurements 

The thermal property of natural abundant and 
perdeuterated HH-P4MeTz were examined by DSC 
(FigJ2J. The thermograms indicate that for both the 
samples there is an endo-thermal peak with a specific 
heat jump near 300K(T C ). This phase transition is due 
to an order-disorder phase transition, namely, partial 
melting of face-to-face 7r-stacking. The similar result 
was observed in an analog of P4RTz 50 . Furthermore, 
another endo-thermal peak appeared at higher tempera- 
ture (>350K) for perdeuterated HH-P4MeTz. The phase 
transition at higher temperature is not due to recrystal- 
lization because an increase of heat of fusion was ob- 
served with an increase of heating rate. It may be due to 
deuteration induced phase transition, but one needs fur- 
ther investigation in order to clarify this phenomenon. 
In this article, we shall give attention only to the phase 
transition near 300K. 

B. 2 H Quadrupolar Echo 

In order to clarify existence of the quasi-ordered phase, 
the 2 H quadrupolar NMR were performed for perdeuter- 
ated HH-P4MeTz at various temperatures (118-393K). 
For the lower temperature phase, in Fig|3| the axially 
symmetric character was observed on the quadrupolar 
powder spectra for the methyl deuteriums. On the other 
hand, the behavior of edge singularities over the transi- 
tion temperature is substantially different from those of- 



ten found in melting of polymers; no averaging (for cases 
of commonly observed melting, very sharp peak at the 
center of spectrum) was observed in the 2 H quadrupolar 
powder pattern and existence of the axially asymmet- 
ric character or another quadrupolar coupling tensor was 
observed. This feature was observed particularly at the 
perpendicular edge of the spectra with a shoulder peak. 
So far, several mechanisms have been pointed out to ex- 
plain the origin of such a slight change of 2 H quadrupolar 
echo spectra. 

First, Hiyama et al. have claimed that the electrostatic 
effect in amino acids causes the quadrupolar coupling 
tensor of deuterium to be asymmetric 74 . It is difficult 
to explain our results by the electrostatic effect, because 
the effect could be more pronounced in the lower tem- 
perature experiments, whereas no asymmetry character 
was observed in the low temperature experiments of HH- 
P4MeTz. Second, Schwartz et al. have pointed out that 
the 2 H- 2 H magnetic dipolar interaction renders the 2 H 
spectrum asymmetric 75 . The effect of dipolar interac- 
tion can be ruled out because of the same reason for the 
electrostatic effect. Third, the asymmetric character is 
also explained by breaking of Cs v symmetry, which have 
been pointed out by Wann et ali^. From Landau's the- 
ory of phase transition, this possibility can be ruled out, 
because a symmetry group for the higher temperature 
phase should belong to a subgroup of a symmetry group 
for the lower temperature phase. Fourth, Kintanar et al. 
have shown that the spectrum can be affected by wagging 
motion and resulted in showing the asymmetrjiZ 7 .. From 
the density measurement and the x-ray diffraction study, 
motion of the methyl groups could be substantially hin- 
dered if exists. However, if partial melting of 7r-stacking 
is present, namely, if a quasi-ordered structure is there, 
the nutation will occur on the methyl group in associa- 
tion with backbone twist. In such a case, the two signals 
derived from the crystalline and the quasi-ordered phases 
can be observed. From the above discussion, it is thought 
that the possibility of the asymmetric quadrupole tensor 
can be ruled out and the two methyl sites with slightly 
different quadrupolar coupling tensor were observed. 



C. Two dimensional Exchange NMR 

2DEX NMR was performed in order to explore more 
detailed molecular dynamics. Figure 0fa) shows the 
2DEX NMR spectra for perdeuterated HH-P4MeTz at 
288K as a function of t m i X . The maximum frequency 
width between cross peaks (Awq) is increased continu- 
ously with the increase of t m ; x (FigQJb)). It was shown 
that this corresponds to the increase of the reorientation 
angle of the methyl group and there exists the motion 
of an order of milliseconds in the crystal. The frequency 
width Acjq is plotted as a function of t m ; x in Fig|5| As- 
suming that there is symmetrical two-site exchange, the 
correlation time(r cx ) can be determined from the fre- 
quency width between cross peaks in the 2DEX NMR 
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spectra by fitting the plot to the following equation 64,66 : 



Awq = Atanh( ), 



(1) 



where A is a fitting parameter. Using this equation the 
correlation time is found with r ex = 6.7ms (t^} = 147Hz) 
at 288K. 

The similar 2DEX NMR spectrum could be also ob- 
tained by the effect of spin-diffusion between deuteriums 
instead of molecular motion^ 8 .. This effect is produced 
by the flip-flop term between the spins which consist of 
magnetic dipole-dipole coupling. In order to estimate the 
spin diffusion time tsd for the zero quantum transition, 
the following equation was considered; 



TSD V27T 



0) 



(2) 



|Wy| is the matrix element of the dipolar flip-flop term 
and a Lorentzian shape is assumed for the resonance 
lines. Indeed, the line width experiments by CPMG 
analog for deuterium shows a Lorentzian shape with a 
FWHM of A/ = 10 Hz. Each term is given by 
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3cos#ij)^ 



(3) 



where is a typical distance between inter-methyl deu- 
teriums and #ij is the angle between the rjj direction and 
magnetic field Ho. The value of | Hij | 2 is estimated as 
560Hz 2 by using a typical distance between the deuteri- 
ums of 0.3nm. 



£ ZQT (w^0) = u m 



(T z Q T rl 



o (T ZQT )- 



hZQT 



(4) 



Here, A/ is assumed to be equal to the homogeneous 
line width(10 Hz). L z( ^ T (ti; — > 0) is estimated as 15ms. 
Eventually, spin diffusion time, tsd, is estimated as 0.3s. 
This value of tsd is about 40~50 times as large as t ox (= 
6.7ms) observed in the measurements. Therefore these 
arguments confirm that the dominant effect which influ- 
ences our 2DEX NMR spectra is not spin diffusion but 
molecular motion. 

The diagonal signals reflect the quantity of immobile 
methyl groups or mobile groups with much faster than 
r mix- The signal intensity on a diagonal, however, di- 
minishes and disappears with increase in t m i X (FigQJc)), 
indicating that no rigid or too rapidly nutating methyl 
group compared with t m j X exists and most of them un- 
dergo very slow motions at 288K. 



D. Proton Longitudinal Relaxation Time 

We performed variable frequency proton spin-lattice 
relaxation time measurements. The proton spin-lattice 



relaxation rate (T^ 1 = Ri) as a function of the frequency 
is shown in FiglHl The relaxation rate shows the depen- 
dence of Ri cx cj -1 / 2 . This relationship between Ri and ui 
is in disagreement with that predicted from the classical 
BPP theory—. The BPP theory in which a fluctuation of 
local field is described as Lorentzian gives the following 
spectral density function, 



J(u) = 



1 + (LUT C y 



(if WT C »1). (5) 



In the region where Ti depends on frequency (ut c 1), 
the relation Ri oc uj~ 2 should be realized. Therefore an 
alternative correlation function should be considered in 
order to explain the dependence obtained by our experi- 
ments. 

So far there are several reports that show the lj^ 1 / 2 
dependence of Ri. The ideas of theories are roughly 
divided into two categories. First, with assuming exis- 
tence of an incommensurate phase, appropriate dynamic 
susceptibility of the classical damped harmonic oscillator 
type for phason branch leads to the lu' 1 / 2 dependence of 
spin-lattice relaxation rat o 82 > 8 ? . At the same time, the 
theory also derives positive proportionality of Ri to tem- 
perature. Then, this theory can be ruled out, because a 
decrease of Ri was observed with an increase of temper- 
ature as shown in Fig[7{a). 

The second theory describes the dynamic susceptibility 
based on a one-dimensional random walk model (details 
of this model are given by the reference). More gener- 
ally, for the weak collision limit {ujt c <C 1), the frequency 
dependence of Ri is different for one-, two-, and three- 
dimensional diffusive motions 80 ' 81 , 



1D:Ri = At 1 J 2 lo- 1 / 2 
2D : Ri = Bt c \u(uj) 
3D : Ri = Ct c - EtI /2 



(6) 



,V2 



The proportionality constants A, B, C and E depend 
on the particular model of hopping motion on a given 
network of atomic sites. It is clear from Fig|S] that one- 
dimensional diffusive motion exists in HH-P4MeTz at our 
measuring frequency. The following equation is obtained 
as a conclusion; 

Ri = M a /(w), (7) 

/H = ^c /2 ^ 1/2 (« < -c 1 ), (8) 

where r c is the correlation time identical to the inverse 
of the diffusion rate and M2 is the second moment of 
the interaction that affects the relaxation. By fitting the 
Eq0 to the result of FigEl the diffusion rate is found to 
be 3±2 GHz. Surprisingly, the rate is an order of 7 larger 
than the exchange rate obtained from the 2DEX NMR 
measurement. 

Since longitudinal relaxation time measurements are 
sensitive to molecular motion with smaller correlation 
time than the 2DEX NMR, it may not be surprising even 



6 



if the two methods (2DEX NMR vs. Ti) show the huge 
difference in the correlation time. 

In the case of undoped trans -polyacetylene, one- 
dimensional electron hopping affects longitudinal relax- 
ation time and diffusion rate is found to be an order of 
10 12 ~I0 13 Hz at room temperature by applying the one- 
dimensional random walk model. It is difhcults in the 
present study to estimate the second moment by taking 
into account electron spin-spin interaction and hyperfine 
contribution, because an amount of free electron per unit 
cell that is proportional to the second moment is not 
known. If mobile electrons govern the Ti measurements, 
the relaxation rate should become large with increase in 
temperature 2 ^. However, as shown in Fig0a), such de- 
pendence was not observed for HH-P4MeTz. Further- 
more, we observed large attenuation of signal intensity at 
diagonal positions in the 2DEX NMR spectra(see above) , 
which implies absense of possible high frequency motions 
sensitive to longitudinal relaxation time measurements. 

Taking into account of absence of high frequency mo- 
tion in the 2DEX measurements and one-dimensional 
fluctuation in proton Ti measurements, we can reach one 
of plausible candidates to explain the huge difference in 
the correlation time, which is wave length dispersion of 
the modulation waves. The distribution of wave number 
produces the localized wave packet by spatial superposi- 
tion. In materials with anomalous dispersion, localized 
waves propagate more rapidly than invdividual waves. It 
can be pronounced that the 2DEX NMR is sensitive to 
the local motion of backbone twist, namely, the nutation 
of individual methyl group was observed (which has a 
character of phase velocity) and that the Ti measuments 
are sensitive to one-dimensional diffusive motion, namely 
the localized modulation waves (group velocity). 

If the relaxation process is governed by thermally acti- 
vated molecular motions with an Arrhenius type, we can 
apply the following expressions; 



E a 



r exp(--) 



(9) 



and if the diffusive motion keeps one-dimensionality in 
our measuring temperature range, from Eqs[7|and|S| 



E 

(i?i) 2 cx t c = r cxp(+-^) 



(10) 



is valid. As Fig[7{b) shows, the correlation time fol- 
lows the Arrhenius law above the phase transition tem- 
perature and the activation energy is found to be 3.4 
kcal/mol. 



increase of echo time(r). In the result of CPMG mea- 
surements in HH-P4MeTz (Fig|S} , we observed the grad- 
ual increase of T 2 " 1 (= R2) for smaller r. The result of 
TPHE measurements was added to the plot since this 
method is thought of as a CPMG method with infi- 
nite echo time (Rj PHE » lim^ m Rf MG ). The decay 
function of TPHE for HH-P4MeTz was exponential with 
first order— and the apparent transverse relaxation time, 
Rj PHE , is much larger than that of the CPMG mea- 
surements. From the arguments by Robertson^, and 
Le Doussal and Sen£&, it is known that when these two 
conditions are realized, the restricted diffusion model can 
be applied; 



i?r HE oc d- 1 



(11) 



For larger r (> 40^s) we also observed the decrease of 
R2 with an increase of r. The behavior of R2 can be un- 
derstood by Sen's theory which treats a restricted diffu- 
sion of magnetization by Bloch-Torrey equation— From 
the theory, the short time regime (r < 40/is) is domi- 
nated by the diffusion process. For the larger echo time, 
the diffusion of the modulation waves reached at bound- 
aries (ex., pinning of modulation waves due to the de- 
fects or interphase between crystal or amorphous and 
the quasi-ordered structure) , and the relaxation behavior 
corresponds to the localization regime, where no diffusive 
character was observed. Therefore the diffusive motion 
in HH-P4MeTz was confirmed by proton transverse re- 
laxation time measurements. Considering the interchain 
face-to-face 7r-stacked structure, we can attribute the dif- 
fusion not to a class of molecular diffusion like chain diffu- 
sion, but to diffusion of the conformational modulation 
waves related to backbone twist. The behavior of the 
temperature dependence of Rj PHE obeys Arrhenius law. 
Therefore the following relation is led from Eallll 



fl™ oc D- 1 



I? - 1 exp(+^), (12) 



and the activation energy was found to be 2.9 kcal/mol 
(Fig|5J). Moreover, ab initio molecular orbital calculation 
with the MP2/6-31G(d) method of TT-4-methylthiazole 
dimer (geometry optimized by the same basis set) shows 
the activation energy of 3.0 kcal/mol for the backbone 
twist (FigHnj). 



13 C CPMAS and 2D Spin-echo NMR 
spectroscopy 



E. Proton Transversal relaxation time 

We report the results of proton transverse relaxation 
time measurements in order to verify whether the diffu- 
sive motion really exists in HH-P4MeTz. One can grad- 
ually incorporate the effect of diffusion into transverse 
relaxation time measured with CPMG method with an 



Figur JTTT a') and (b) show the 13 C cross-polarization 
and magic-angle sample spinning(CPMAS) spectra for 
HH-P4MeTz at 293K. For the spectrum, only signals de- 
rived from methyl carbons appear over the shift region 
of 15-21 ppm. It is worth to point out that the shoul- 
der signal of the methyl carbon are observed at 17 ppm. 
Although the previous X-ray diffraction study^ suggests 
that there exists an electronically unique methyl carbon 
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in the unit cell of HH-P4McTz, the NMR spectrum indi- 
cates that there exist at least two electronically distinct 
methyl carbons. 

Chemical shielding tensor(CST) is more informative 
than its isotropic shielding, about three-dimensional elec- 
tronic structure. It is difficult, unfortunately, to deter- 
mine principal values of CST for nuclei with small CSAs, 
such as methyl carbons. In order to gain further in- 
sights on the chemical shielding difference, we performed 
solid-state 13 C two-dimensional spin-echo(2DSE) CP- 
MAS measurements and extract the information about 
principal values of CST for the methyl carbons. Fig lllf c) 
and (d) shows the solid-state 2DSE spectra and the 
best fitted simulated spectra for the methyl carbons in 
HH-P4MeTz. The values of error parameter, e, were 
determined as 8.3xlCT 3 for the Mei of HH-P4MeTz, 
3.4xl0~ 2 for the Me 2 of HH-P4MeTz, respectively. Ta- 
ble 0] summarized the observed principal components 
of chemical shielding tensor for the methyl carbons in 
HH-P4MeTz. It is realized that the chemical shielding 
anisotropies for the two distinctive signals of the methyl 
carbons in HH-P4MeTz are quite different from each 
other. Each principal components for the Mei signal 
is less shielded compared to those for the Me2 signal. 
This can be mostly pronounced for the midfield compo- 
nent, S22', the difference is approximately 6ppm. This 
dominates the other principal components and gives the 
isotropic shielding difference of 67%. We found that the 
inhomogeneously broadened signals of the methyl car- 
bon is due to rearrangements of the backbone torsion 
in the crystalline lattice, namely, due to coexistence of 
the quasi-ordered phase with the crystalline phase, in 
which the polymer chains stack in a manner similar to the 
crystalline polymers but the tt-tt interactions are weaker. 
With the consideration of the results of 2DEX, the quasi- 
ordered phase might dynamically coexist with the crys- 
talline phase. 

The experimental results of 13 C CST is further under- 
stood by ab initio chemical shielding calculation. In Ta- 
bleman averaging of the shielding tensor with respect to 
P was carried out under the consideration of the Boltz- 
mann factor. The experimental chemical shielding for 
the methyl carbon in HH-P4McTz is qualitatively repro- 



duced by variation of backbone twist in the quasi-ordered 
phase. 



V. CONCLUSION 

Solid-state dynamics of regio-regulated HH-P4MeTz 
was investigated by various solid-state NMR spec- 
troscopy. DSC measurements show that both natu- 
ral abundant and perdeuterated HH-P4MeTz causes the 
phase transition near 300K. In the 2 D quadrupolar NMR 
spectra at various temperatures also suggest that the 
phase transition from the crystalline to the quasi-ordered 
phase occurs. The two-dimensional exchange NMR ex- 
periments revealed that the motion of methyl groups with 
an order of milliseconds exists in HH-P4MeTz and the ex- 
change rate of such the motion was estimated as r^ 1 = 
147 Hz. Furthermore, the relationship of Ri ~ lu~ x / 2 was 
observed from proton longitudinal relaxation time mea- 
surements, which might be due to the one-dimensional 
diffusion-like motion of conformational modulation waves 
related to the nutation of methyl groups along the chain. 
Existence of diffusion of modulation waves was also con- 
firmed from proton transverse relaxation time experi- 
ments. We obtained the diffusion rate of 3±2 GHz by 
calculating the dynamic susceptibility with assumption of 
the ID random walk model. The discrepancy between the 
results of the Ti measurements and 2DEX NMR may be 
due to the anomalous dispersion of the modulation wave. 
13 C CPMAS NMR measurements, an analysis of princi- 
pal components of CST for the methyl carbons deter- 
mined by the 2DSE measurements, and chemical shield- 
ing calculation show that there exists a quasi-ordered 
phase in HH-P4MeTz and that the backbone twist is 
highly correlated with the methyl nutation. 

From the 2 D 2DEX NMR and the *H Ti measure- 
ments, it is also found that HH-P4MeTz has an anoma- 
lous dispersion of the modulation wave. This research 
is the first example that have measured phonon disper- 
sion by NMR. Such dynamics may be detectable due to 
existence of distinctive one-dimensional fluctuation in tt- 
conjugated systems. 
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TABLE I: Principal components of observed nuclear shielding 
tensors of the methyl carbons in HH-P4MeTz (units in parts 
per million from tetramethylsilane.) 



material 


^iso 


Su 


822 S33 


8 


V 


n a 




e 


HH-P4MeTz(Mei) 


20.1 


29 


23 8 


-12 


0.5 


21 


0.41 


8.3xl0~ a 


HH-P4MeTz(Me 2 ) 


17 


28 


17 6 


-11 


1.0 


22 


0.00 


3.4xl0~ 2 



a) span: 57 = 8u — S33 
b) skew: k = 3(^22 — &so)/(5ii — 833) 

TABLE II: Principal components of calculated absolute nu- 
clear shielding tensors of methyl carbons in HH-P4MeTz 
(units in parts per million.) 



model 


Method 


Tiso 


(Til 


U22 a 33 


8 r] 


Q 


K 


s- trans 


SCF 


186.7 


178.2 


180.4 201.4 


-14.7 0.150 


23.2 


0.8 14 


averaged 


SCF 


189.8 


180.3 


184.5 204.4 


-14.6 0.288 


24.1 


0.6 60 


s- trans 


MP2 


190.6 


181.2 


184.4 206.3 


-15.7 0.204 


25.1 


0.7 41 


averaged 


MP2 


194.2 


184.1 


188.7 209.8 


-15.6 0.295 


25.7 0.6 42 
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FIG. 1: Chemical structure and schematic representation of 
the alternative layered molecular packing in HH-P4MeTz; 
Neighboring molecules take a face-to-face packing (b axis). 
The X-ray diffraction study shows that thiazole rings recur- 
ring along c axis take coplanar structures and form a highly 
extended 7r-conjugated network. It takes the alternative lay- 
ered structure: a layer constituted by highly densed methyl 
groups with two dimensional manner and a layer with the 
face-to-face 7r-stacking. 



d 3 -HH-P4MeTz, 5K/min. 
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FIG. 2: DSC charts for powder samples of perdeuterated HH- 
P4MeTz and natural abundant HH-P4MeTz. The heating 
rate is shown in the figure. After the first heating scan, the 
samples were quenched by liquid nitrogen. 
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FIG. 3: 41MHz 2 D duadrupolar echo NMR spectra for the 
powdered sample of perdeuterated HH-P4MeTz at various 
temperatures (118 ~ 393K). The dashed line is only a guide 
for eyes. The axially symmetric spectra for the methyl deu- 
terium with C3t, rotation were observed below 300K. Over 
300K, the shoulder peaks were detected near the parpendic- 
ular edge singularities. 
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FIG. 4: The two-dimensional exchange 2 D NMR spectra for 
HH-P4MeTz at 288K with various mixing time (a), The cross- 
diagonal (b) and the diagonal (c) slice of the 2DEX spectra for 
HH-P4MeTz. Tightly fixed or rapidly moving methyl groups 
will maintain intensities of the diagonal peaks if exist. In 
our measuring time scale, such the maintainance was not ob- 
served, indicating the most methyl groups undergo the slow 
motion associated with backbone twist. 
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FIG. 5: The variation of the maximum frequency difference, 
Aujq, with mixing time (t m i x ) for the non-diagonal peaks. 
The dotted line is a theoretical curve of Eq0 
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FIG. 6: The frequency dependence of proton longitudinal re- 
laxation rate for HH-P4MeTz at 288K. Each symbol shows 
the results of three different measurement techniques: the 
saturation-recovery method (closed circle), the spin- locking 
method (closed triangle) and the Jeener-Broekaert method 
(closed square) respectively. The dotted line is the relation 
of Ri <x a;" 1 / 2 . The wide range of frequency of the rotating 
frame was carried out by using a probe with a micro coil. 
The frequency of the dipolar order was determined by the 
Van Vleck's second moment method. 
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FIG. 7: The temperature dependence of proton relaxation 
rate over the temperature range from 290 to 355K. The value 
of Ri becomes smaller with elevation of temperature (a). The 
plot of Ri as a function of 1000/T (b). When one- dimensional 
diffusion is responsible for the behavior of Ri, R^ is propor- 
tional to t c (see Eg 1101 . The thermally activated process with 
an Arrhenius type was able to be observed. 




Echo Time 



FIG. 8: The dependence of transverse relaxation rate on 
CPMG echo time. The data for two-pulse Hahn echo(TPHE) 
were added to the plots since the TPHE method is thought as 
a CPMG experiment with infinite echo time. The solid lines 
denotes the extrapolation of the transverse relaxation rate at 
the short time regime to the relaxation rate of TPHE. 
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E a = 2.9 kcal/mol 
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FIG. 9: The variation of R2 by TPHE as a function of temper- 
ature near above T c . The thermally activated process was ob- 
served and the activation energy was found to be 2.9 kcal/mol. 
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FIG. 10: A model compounds for quantum mechanical chem- 
ical shielding calculation (a). Chemical shieldings for the 
carbons marked with an asterisk were calculated. The two 
methyl carbons in TT-P4MeTz is electronically equivalent. 
The MP2 energy for HH-P4MeTz as a function of backbone 
twist of (3 (b). 
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FIG. 11: 67.9MHz 13 C cross-polarization and magic- 
angle sample spinning(CPMAS) spectra for HH-P4MeTz at 
293K(a,b). (c),(d) 100.0 MHz observed and best fitted simu- 
lated solid-state 13 C two dimensional spin-echo CPMAS spec- 
tra for the methyl carbons in HH-P4MeTz(Mei ;c) and HH- 
P4MeTz(Me 2 ;d). 



